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Is there a compact companion orbiting the late O-type binary star 
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ABSTRACT 

We present a multi-wavelength (X-ray, y-ray, optical and radio) study of HD 164816, a late 
O-type X-ray detected spectroscopic binary. X-ray spectra are analyzed and the X-ray pho- 
ton arrival times are checked for pulsation. In addition, newly obtained optical spectroscopic 
monitoring data on HD 164816 are presented. They are complemented by available radio 
data from several large scale surveys as well as the FERMI y-ray data from its Large Area 
Telescope. We report the detection of a low energy excess in the X-ray spectrum that can be 
described by a simple absorbed blackbody model with a temperature of ~ 50 eV as well as a 
9.78 s pulsation of the X-ray source. The soft X-ray excess, the X-ray pulsation, and the kine- 
matical age would all be consistent with a compact object like a neutron star as companion 
to HD 164816. The size of the soft X-ray excess emitting area is consistent with a circular 
region with a radius of about 7 km, typical for neutron stars, while the emission measure of 
the remaining harder emission is typical for late O-type single or binary stars. If HD 164816 
includes a neutron star born in a supernova, this supernova should have been very recent and 
should have given the system a kick, which is consistent with the observation that the star 
HD 164816 has a significantly different radial velocity than the cluster mean. In addition we 
confirm the binarity of HD 164816 itself by obtaining an orbital period of 3.82 d, projected 
masses raj sin 3 / = 2.355(69) M , m2sin 3 / = 2.103(62) M apparently seen at low inclination 
angle, determined from high-resolution optical spectra. 

Key words: stars: neutron, stars: individual: HD 164816, stars: individual: 
2XMM J180356. 8-241845. 



1 INTRODUCTION 

The expected total number of Neutron stars (NSs) in our Galaxy is 
predicted to be 10 8 -10 9 (Narayan & Ostriker 1990) of which iso- 
lated NS should form the majority. Until today there are ~2000 
known radio pulsars and only seven known isolated thermally emit- 
ting NS (INS) called the Magnificent Seven (Haberl 2007, Kaplan 
et al. 201 1). Since the discovery of the first INS (RX J1856.5-3754) 
in 1996 (Walter et al. 1996) the search for more thermally emit- 
ting NSs is an ongoing process. Those seven objects have been 
recognized by their high X-ray to optical flux ratio and their rather 
soft X-ray emission represented by low X-ray hardness ratios. Thus 
looking for objects with similar properties is one way to find new 
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candidates, see e.g. Pires et al. (2009). It is clear that many candi- 
dates can be missed, for example when they are harbored in binary 
or multiple star systems, since the X-ray flux is dominated by the 
host star, or as compact companions to runaway stars (Posselt et al. 
2008). Searches for Pulsar companions around OB runaway stars 
have been performed by Philp et al. 1996, Sayer et al. 1996. 
We search for X-ray pulsations from all X-ray sources near or 
identified with galactic OB stars, in search of non-interacting (and 
hence effectively isolated) NSs that remain bound in a stellar sys- 
tem following the supernova. About 10% of such systems remain 
bound after the first supernova (Kuranov et al. 2009). 
HD 164816 is located in the direction of the young open cluster 
NGC 6530. Prisinzano et al. (2005) recently investigated in depth 
the distance of this cluster: Literature values range from 560 pc 
(Loktin & Beshenov 2001) to 2000 pc (Walker 1957, van den An- 
cker et al. 1997). The mean of the distance estimates (in table 1 
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of Prisinzano et al. (2005) and the value found by Prisinzano et al. 
(2005) themselves, 1250 pc) is 1543 ± 345 pc, which we now use 
in this paper; all but one (560 to 71 1 pc, Loktin & Beshenov 2001) 
values found in the literature are consistent with this value within 
1 cr. Within 1 cr, our value is also consistent with the distance of 
~ 1250 pc obtained by Prisinzano et al. 2005. 

HD 164816 itself is located at a = 18h 03m 56.866s and 
6 = -24°18'45.22" and has V=7.09 mag (Reed 2003). Recently, 
Megier et al. (2009) determined the distance towards HD 164816 
by interstellar Ca II absorption to be 864 ±136 pc. Hence, the dis- 
tances of the NGC 6530 cluster and the star HD 164816 are deviant 
by 1 to 2 cr, which may indicate that HD 164816 is somewhat fore- 
ground to the cluster, to be discussed below. 
The region around it is also known as the Lagoon Nebula (NGC 
6523 = M 8) which is one of the brightest HII regions in the Galaxy 
where star formation has started a few 10 7 yrs ago (van den Ancker 
1997). However, the membership of HD 164816 is uncertain. 
NGC 6530 has about 2000 known members which display an age 
gradient (Damiani, Prisinzano, Micela & Sciortina 2006). From 
March 2001 to July 2003 HD 164816 has been in the field of view 
of seven X-ray observations, one with XMM Newton and six with 
the Chandra X-ray observatory (CXC). 

In the following we will present the X-ray spectral and temporal 
properties from XMM PN and Chandra ACIS-I data indicating a 
compact companion candidate (§2). We present the orbital parame- 
ters of the spectroscopic binary O-type star in §3. Then we present 
the available radio survey data as well as the analysis of y-ray data 
obtained by FERMI 's Large Area Telescope (EAT) in §4. A discus- 
sion of the physical nature of the possible companion will be given 
in §5 and concluding remarks will be made in §6. 



2 X-RAY DATA ANALYSIS 

We report here on the spectral and temporal analysis of the X-ray 
source 2XMM J180356.8-241845 that is coinciding with the opti- 
cal position of HD 164816 and has been found by our search for 
pulses in X-ray sources near galactic O-type stars. 
We extracted events from the archived pipeline produced XMM PN 
data (see table [TJ in a circular region with a radius of 30 arcsec 
around the centroid of 2XMM J180356.8-241845 located approxi- 
mately 3. 1 arcmin from the pointing of the observation by using the 
XMM Science Analysis System (SAS) version 10. As background 
we extracted events in a nearby source free region of similar size. 
We used only single-pixel event types for the spectral analysis as 
well as for the performed source detection in the first energy band 
of XMM-Newton (0.2 - 0.5 keV) as there are known calibration is- 
sues. Event types higher than single-type are known to result in an 
elevated soft X-ray background which is most prominent in the re- 
ported first energy bancQ 

Extracting the background region is a somewhat complicated task 
as the field is very crowded and contaminated by several fringes. 
Although the thick filter was in place in the XMM Observation, 
those might result from incomplete blocking of UV/optical light 
coming from the 5th magnitude star 9 Sgr (see Fig[T|( or possibly 
originating from stray light of the nearby Low Mass X-ray binary 
(LMXB) GX 9+1 (Langmeier, Sztajno, Trumper & Hasinger 1985) 



around 1 deg outside the field of view (fov). Checking the back- 
ground light curve of this observation we could not identify any 
times of elevated background and thus can use the whole observa- 
tion time as a good time interval. Nevertheless we also investigated 
the spatial distribution of the arriving photons in the background 
extraction region and could not find any significant concentration, 
rather a homogeneous distribution (see inset in Fig[T]left panel). 
The Chandra observations do not show any fringes or an otherwise 
enhanced background as the threshold for optical contamination is 
V ~ 3.1 for the used ACIS-I detecto^](HD 164816 has V=7.09). 
Hence, the chosen background regions are not contaminated by any 
source or fringe. 

2.1 Spectral Analysis with XMM Newton 

As a first step we have re-analyzed the X-ray spectrum of 
2XMM J180356.8-241845. Its spectral properties have already 
been reported in Rauw et al. 2002 where they performed an X-ray 
population study on an XMM Newton observation and among other 
things conclude that the X-ray emission of HD 164816 is typical 
for an 09.5 III-IV star. However, they did not take the low energy 
part of the spectrum from 0.2 - 0.5 keV into consideration. 
Trying to reproduce their results we first re-analyzed the XMM PN 
dataset of ObsID 0008820101 carried out with the thick filter in use 
and centered at 9 Sgr. 

As a following step Response and Anxillary files have been com- 
puted by using the XMMSAS tasks RMFGEN and ARFGEN. The 
data then have been grouped by a minimum of 20 counts per en- 
ergy bin in order to minimize the scattering of data points. After 
background subtraction there are 961 net source counts available. 
Including as well the first XMM Newton energy band (0.2 - 0.5 
keV) in the spectral fit by using a warm absorbed MEKAL plasma 
(Mewe et al. 1985, Kastra 1992) results in an unsatisfying fit of x' v 
equal 1.57 with 34 d.o.f which is not in agreement with what Rauw 
et al. 2002 have found (i.e. ^=1.08 in 0.5 - 5 keV). By looking 
at the spectrum we clearly notice an excess in the range 0.2 - 0.5 
keV that cannot be modeled by the used absorbed MEKAL plasma 
model. In addition we recognize some feature around 0.3 keV but 
we ignored data at the energy range 0.27 to 0.32 keV as there is 
a drastic dip in the effective area of the EPIC PN detector in that 
energy range. 

We have carried out a source detection in the energy range 0.2- 
0.5 keV in order to exclude an artificial nature of the detected low 
energy excess. This energy part is believed to be mostly from the 
possible compact companion and has not been regarded by Rauw 
et al. (2002). We detected 2XMM J180356.8-241845 in this en- 
ergy range with a detection maximum likelihood of L ~235 with 
L = - In p where p is equal to the probability that the detected 
signal was generated by a random fluctuation. Since we only took 
single-pixel events into consideration we can rule out that the ex- 
cess is due to elevated background. As the given likelihood value 
for the source detection can be considered significant we conclude 
that the excess is of real nature (see right panel in Fig[TJ. A source 
detection in the energy range 2.0 to 10 keV leads to a null detec- 
tion, hence the source is soft. Improving the fit by blindly adding 
a blackbody component could not be achieved thus we decided to 
use a more sophisticated approach. Therefore we tried to subtract 
the MEKAL component from the spectrum. 



1 Please see document XMM-SOC-CAL-TN-0068 
http://xmm2.esac.esa.int/extemal/xmm_sw_cal/calib/documenta- 
tion/index.shtml for reference 
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Figure 1. Left Panel: Image showing a 17x 14 arcmin 2 close-up view of the EPIC PN observation centered at the position of HD 164794 (9 Sgr) in the energy 
range 0.2 - 10.0 keV. Shown are the position of HD 164816 with the circle indicating the 30 arcsec extraction region together with the background extraction 
region. As well seen is the rather large contamination by other sources in the field and the contribution of fringes from the optical emission of 9 Sgr. In addition 
surface plots of the extracted source (left inset) and background regions (right inset) are shown. Right Panel: Here we show HD 164816 in the energy range 
0.2-0.5 keV with only single events displayed. 



Table 1. X-ray observation log of the used XMM-Newton and Chandra observations 



Observation ID 


Pointing 


Start Time 


exposure time (ks) 


Detector 


offset from HD 164816 (arcmin) 




XMM-Newton 










0008820101 


HD 164794 (9 Sgr) 


2001-03-08 11:21:27 


23.572 


MOSl/2,PN 


3.07 




Chandra 










3754 


M8 


2003-07-25 17:27:12 


129.600 


ACIS-I 


4.27 


4397 


M8 


2003-07-24 10:07:26 


14.820 


ACIS-I 


4.27 


4444 


M8 


2003-07-28 00:00:07 


30.190 


ACIS-I 


4.27 



In order to eliminate the thermal plasma component we used 
XSPEC (version 12.6) to simulate an absorbed MEKAL plasma in 
the energy range 0.5-2.0 keV as there is negligible information in 
the range 2.0-10.0 keV. As an input the best fit parameters for tem- 
perature and absorption reported in Rauw et al. 2002 were used. 
Afterwards this simulated spectrum is used to create a fake data 
set with the same exposure time, background, response and anxil- 
lary files as used for the observation of HD 164816. The faked data 
are grouped according to the EPIC PN data (i.e. a minimum of 20 



counts/bin). This grouped MEKAL spectrum is then used as back- 
ground for the observation in order to subtract the contribution of 
the thin plasma model. 

Visually checking the produced spectrum confirms the MEKAL 
component to be consistent with zero. The remaining excess is fit- 
ted with an absorbed blackbody model which leads to a best fit 
temperature of around 53.7 eV that now serves as an initial input 
for the future complete analysis. As a next step the entire spectrum 
is fitted with an absorbed thin plasma model plus the blackbody 
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HD1648'6 with abs(BB + abs MEKAL) model Coadded CXC ACIS-I spectrum for HD1648'6 




Figure 2. Left Panel: Shown here is the best fit absorbed blackbody plus MEKAL model (continuous line) to the XMM PN data together with the residuals 
for a blackbody temperature of kTy, ~48.79 eV. Data points from 0.27 keV to 0.32 keV are not taken into consideration in the fit due to the drastic dip in 
the effective area around 0.3 keV. The dash-dotted line represents the best fit absorbed MEKAL model to the data and one can see the already mentioned 
excess from 0.2 to 0.5 keV. Right Panel: The coadded Chandra ACIS-I observations are shown here together with the best fit absorbed blackbody + MEKAL 
plasma model and their corresponding residuals for a blackbody temperature of kT^j, ~49.85 eV. Likewise the dash-dotted line represents the best fit absorbed 
MEKAL model, where the excess from ~0.3 to ~0.5 keV becomes visible. 



component by iteratively fixing and freeing the parameter pairs for 
the blackbody and MEKAL plasma model (i.e. pairs of temperature 
and normalization) which lead us to a best fit blackbody tempera- 
ture of 52.36+™! eV with x l = 1-06. 

As the interstellar column density derived by Diplas & Savage 
1994, 0. 15+q^-IO 22 cirr 2 , is significantly smaller than the value 
derived by our fit (0.40+qqj-IO 22 cm -2 ) we decided to introduce a 
second systemic absorption component. Such a systemic absorp- 
tion was already suggested by Rauw et al. 2002 and could pos- 
sibly be attributed to a stellar wind. By introducing such a com- 
ponent we can put tighter constraints on the errors of the individ- 
ual parameters and get an interstellar absorption of 0.08^ JJj-lO 22 
cnT 2 which is consistent with the one found by Diplas & Sav- 
age 1994. The resulting blackbody temperature is 48.79+^5 eV 
(see Fig[2]l. Values for the systemic absorption and MEKAL tem- 
perature are 10 22 cirr 2 and 0.24±0.02 keV respectively. 
The absorbed flux in the energy range 0.2-2.0 keV of 1.08xl0~ 13 
erg cnr 2 s~' is found to be higher than the presented value for 
the energy range 0.5-5.0 keV by Rauw et al. 2002. For complete- 
ness we have computed as well the unabsorbed flux in the 0.2-2.0 
keV range (2.89xl0~ 13 erg cnr 2 s~'). In addition we have used the 
blackbody model normalization to give an estimate on the radius 
of the emitting area. Using a distance of 864±136 pc results in a 
radius of the (circular) emitting area of 7. 14± 1. 12 km which is in 
the range of a NS (see Hohle et al. 2012). Adopting the same dis- 
tance range we used the MEKAL normalization to give an estimate 
on the emission measure (EM) of HD 164816, which turns out to 
be 3.93+0. 10xl0 54 cirr 3 . Those values are typical for 09V single 
stars (see e.g. Bhatt et al. 2010 and Naze et al. 201 1). An overview 
of the most important parameters can be found in table[2] 



2.2 Spectral Analysis with Chandra 

HD 164816 has been in the field of view of six Chandra ob- 
servations of which we unfortunately had to exclude two as 
the high energy transmission grating (HETG) was in place and 
HD 164816 thus not covered by the active part of the detector. 
The remaining four observations (ObsIDs 977, 3754, 4397, 4444; 



see table [TJ taken by the Advanced CCD Imaging Spectrometer 
in Imaging mode (ACIS-I) have been reprocessed individually, 
as follows. First we extracted events in a circular region around 
HD 164816 with a radius of 10 arcsec for each observation and 
as there are neither any fringes nor any other signs for an elevated 
background present we arbitrarily chose one background field per 
observation with a radius of 20 arcsec in a source free region in 
close proximity to HD 164816. We chose a somewhat bigger ex- 
traction region for the background as the background level of Chan- 
dra observations is in general very low. 

Response and Anxilary files have been created afterwards for 
source and background regions each by using the CIAO (version 
4.2) tools MKRMF and MKARF. As we noticed that the source 
counts are low in the single observations we inspected the effec- 
tive area of each observation and found the differences for ObsIDs 
3754, 4397 and 4444 to be indistinguishable and hence combined 
these three event files. However for ObsID 977 we find the effec- 
tive area to be significantly smaller and thus did not include this 
observation. Coadding the remaining three observations lead to a 
total effective exposure time of ~ 172 ks with 1203 counts available 
for HD 164816 and thus having slightly higher statistics compared 
to the XMM observation. 

For the following spectral analysis we applied the same methods 
and constraints as to the XMM PN observation but we used the 
C-statistic instead of x 1 m that case. As a best fit blackbody tem- 
perature we get 49.85+j 2 !? eV with C = 168.77 for 109 d.o.f. in the 
energy range of 0.3-2.0 keV (see Fig(2] right panel), as the ACIS- 
I detector is only calibrated from 0.3 to 11.0 keV. Comparing the 
MEKAL plasma temperature and both absorption components ob- 
tained from the Chandra data with the values from XMM PN data 
we note that they are in agreement. This is also true for the values 
of the MEKAL plasma temperature in the Rauw et al. 2002 anal- 
ysis. The computed absorbed and unabsorbed flux are 1.19xl0~ 13 
erg cnr 2 s~' and 2.19xl0~ 13 erg cnr 2 s~' respectively. Computing 
the radius of the emitting area 864±136 pc we end up with K emi , 
= 7.07± 1.11 km which is consistent with the value derived from 
the XMM fit. By using the MEKAL normalization factor we find 
the EM in the range 3.76±0.09xl0 54 cnT 3 , which is similar to the 
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Table 2. Spectral parameters inferred from fitting the Chandra and XMM-Newton spectra of HD 164816. 



Model " 


X 2 1 C-stat 


D.O.F. 


interstellar nn / systemic tin 
(10 22 cm" 2 ) 


kl MEKAL 

(eV) 


kTsB 
(eV) 


R 

(km) 


fx (0.2-2.0 keV) 
(10~ 13 erg cm" 2 s _1 ) 








XMM-Newton 


(0.2 - 2.0 keV) 








abs*MEKAL 


X 2 = 53.30 


34 


35 +0 08 

U - J -0.07 


237.81+j^' 






7.90 


abs(BB+ abs*MEKAL) 


X 2 = 32.84 


31 


0.08^/0.32^ 


237.81^ 


48.79^ 


7.14±1.12 


2.89 


Model " 


X 2 1 C-stat 


D.O.F. 


interstellar nn / systemic nn 
(10 22 cm 4 ) 


kT MEKAL 

(eV) 


kT BB 
(eV) 


Remit 

(km) 


fx (0.3-2.0 keV) 
(10~ 13 erg cm -2 s ) 








Chandra 


(0.3 - 2.0 keV) 









abs*MEKAL 184.68 112 °- 39 -002 227.70+^° — — 1.01 

abs(BB+ abs*MEKAL) 168.77 109 0.09^^/0.28+°^ 238.90+[ 2 ^ 2 49.85+^ 2 7.07±1.11 2.19 



BB = blackbody; MEKAL = hot diffuse gas model based on Mewe et al. 1985, 1986, Kaastra 1992 and Liedahl et al. 1995 



XMM data. 

We summarize our X-ray results in Table[2] 



2.3 X-ray pulsation search with XMM Newton 

We performed a periodicity search at the position of 
2XMM J180356.8-241845. For performing the pulsation search 
only the XMM Newton observation has been reprocessed as the 
EPIC PN detector was used which is the only available instrument 
for this source with a sufficiently short readout time (73 ms in 
that configuration). In particular there are 1284 arrival times in 
this set that have been searched for a periodic signal by means 
of the bayesian approach developed by Gregory & Loredo 1992 
for the detection of a signal with unknown shape. Applying this 
method leads to a detection of a sie nal at 9.7804+°;°*™ s (see Fig[3} 
with a probability that the data favor a periodic model over a 
constant model to be 0.56. This resulting period could afterwards 
be reproduced by applying the Z 2 test with m being the number 
of harmonics (Buccheri et al. 1983) in the period range 1 to 10s 
with a stepsize equal to the independent fourier spacing (IFS) i.e. 
l/°T spm , This number of harmonics has been optimized by using 
the H-Test (De Jager, Swanepoel & Raubenheimer 1989) and was 
found to be maximized for m=l. The resulting Z 2 value is equal 
to 22.67 at a period of 9.7813+q5Joo7 s w ' m me number of expected 
peaks exceeding that Z value to be 5.25x 10~ 6 . The given errors are 
1 cr errors. Finally a pulsation search with a method based on the 
cash statistic (Cash 1979) as described in Zane et al. (2002) has 
been applied to the arrival times and could determine the period to 
be 9.78±0.06 s (again in the range 1 to 10s and using the IFS). The 
errors are again at the 1 cr level. 

Hence, the period found by the bayesian method, being the proper 
choice for this kind of data as it uses an unbinned approach, is 
confirmed by the cash and Z 2 -test. In addition we have created a 
simulated data set of equal length in exposure time with 1284 ran- 
domly distributed photons. Applying the Z 2 -Test and the cash-test 
to that set of data lead to no detection at ~9.78 s. Furthermore the 
most likely periods found by those two tests in the faked data set 
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Figure 3. Here we show the periodogram of the detected ~9.78 s pulsa- 
tion for the X-ray source identified with HD 164816 found by the bayesian 
approach with XMM Newton. 



are not in accordance (P z 2 ~5. 15 s and Y cash ~ 1 . 12 s). 
The CXC ACIS-I data sets although improving the number of 
source photons have not been used for timing analysis as the ACIS 
frame time is ~3.24 s and thus too long to significantly detect such 
a period. 

The phase-folded X-ray light curve for 2XMM J180356.8- 
241845 in the energy range 0.2-2.0 keV has a pulsed fraction 
of 60.4±15.4 % (c.f. Fig|4](. As a final step we extracted arrival 
times for the brightest X-ray source in the field of view (i.e. the 
04 star 9 Sgr) in the energy range 0.2-2.0 keV which results in 
18769 counts. We folded the created light curve with the obtained 
period which results in a constant non-periodic signal. This is 
likewise true for the extracted background region, again in the 
energy range 0.2-2.0 keV. Considering this non-periodic signal of 
the background and of 9 Sgr, we conclude that the observed period 
in HD 164816 is not due to some background modulations caused 
for example by the readout time or multiples of it. 
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HD 164816 in 0.2 - 2.0 keV into 9.78s binning=0.1 




0.5 1.0 1.5 

Phose 



Figure 4. In this figure we show the phase-folded light curve for the X-ray 
source in the energy range 0.2-2.0 keV folded into the resulting best period 
of 9.78 s resulting from XMM Newton data. For clarification we show two 
full cycles with a binning of 0. 1 . 



Both the soft excess and the ~10 s pulsation would be typical for 
a NS in the system assuming that it will not be interacting with 
its host star and can hence be treated as effectively isolated (see 
Haberl 2007 for typical periods of thermally emitting isolated NS), 
namely ~3 to ~ 1 1 s. 



3 ORBITAL PARAMETERS OF HD 164816 
3.1 Optical data analysis 

Optical echelle spectroscopy of HD 164816 has been obtained 
at two observatories in Chile: the Cerro Armazones Observatory 
(OCA, about 20 km east of the Paranal Observatory) and at the 
Cerro Tololo Interamerican Observatory (CTIO). 
Optical spectroscopy at the OCA was secured using the Bochum 
Echelle Spectrographic Observer (BESO) fiber-fed from the focus 
of the 1.5m Hexapod Telescope (see Fuhrmann et al. 2011). Fif- 
teen spectra were obtained between 2009-04-21 and 2010-03-28 
covering the wavelength range of 3530 - 8860 A with a spectral 
resolution of R = 50000 (Fuhrmann et al. 201 1) (see Table[5]l. The 
data were reduced using dedicated scripts written under MIDAS. 
The reduction includes overscan, bias and flatfield correction; in- 
dividual echelle orders were extracted, wavelength calibrated and 
normalized to the continuum. Finally, cosmic spikes have been re- 
moved. The cleaning was complicated by the fact that most He 
and H lines contain narrow (and variable) emission features. More- 
over, for each night there is just one exposure. Wavelength cali- 
bration has been improved using telluric bands close to 6900 and 
7600 A and the spectrum of rj CMa as a template. The improved 
radial-velocity system is stable to about 100 m/s. The zero point of 
the radial-velocity system was checked by measuring the IAU stan- 
dard star HIP 910, which is a slow rotator, in 6 nights. The average 
radial-velocity is 14. 16+0. 13 km/s which is in good agreement with 
14.4+0.9 km/s found by Evans (1967). 

Optical spectra at the CTIO were obtained using the high disper- 
sion optical Blanco echelle spectrograph fiber-fed from the 1.5m 
SMARTS telescope. Fourteen spectra taken between 2010-07-25 
and 2010-09-14 cover the range from 4820 - 7120 A (see Table 
[3]>. The R «20000 spectra were extracted using software written 



Table 3. Journal of spectroscopic observations of HD 164816 obtained at 
the CTIO and OCA observatories 



Spectrum 


HJD 


Exp. 


Spectrum 


HJD 


Hxp. 




2 400 000+ 


[sec] 




2 400000+ 


[sec] 


OCA.20090414 


54936.8391 


1800 


OCA_20 100328 


55284.8365 


1500 


OCA.20090421 


54943.8225 


1800 


CTIO.20100725 


55403.7934 


200 


OCA.20090508 


54960.8685 


1800 


CTIO.20100728 


55406.7316 


300 


OCA_20090707 


55020.8242 


1800 


CTIO_20100731 


55409.6988 


300 


OCA.2009I010 


55115.5480 


1800 


CTIO.20100804 


55413.6327 


300 


OCA.2009I01S 


55123.5011 


1800 


CTIO_20100805a 


55414.6293 


300 


OCA.2009I020 


55125.5094 


1800 


CTIO.20100805b 


55414.7672 


300 


OCA_2009I024 


55129.5112 


1800 


CTIO_20100812a 


55421.6210 


600 


OCA_2009I025 


55130.5172 


1800 


CTKX20100812b 


55421.7031 


600 


OCA.20091027 


55132.5113 


1800 


CTIO.20100820 


55429.6515 


400 


OCA.20100321 


55277.9383 


1500 


CTIO.20 100904 


55444.5716 


600 


OCA.20100323 


55279.9109 


1500 


CTIO.20 100907 


55447.5156 


600 


OCA_20100324 


55280.8909 


1500 


CTIO-20100908 


55448.5887 


600 


OCAJ0100325 


55281.8479 


1500 


CTIO_20 100909 


55449.5630 


600 


OCA.20100327 


55283.8057 


1500 


CTIO.20100914 


55454.5746 


600 



in IDL|j The orders are traced using the cross-dispersed flat im- 
ages. The spectrum is extracted using a boxcar extraction, as is the 
flat spectrum. The spectrum is then divided by the extracted flat 
spectra. In general we obtained 3 spectra at each epoch. These are 
scaled, median-filtered to reject cosmic rays, and summed. Wave- 
length calibration is based on a Th-Ar calibration lamp exposure 
taken just prior to each stellar observation. The zero-point of the 
wavelength scale is uncalibrated at the level of 1 pixel (3 km/s). No 
attempt has been made to convert to a flux intensity scale. 

Although the zero point of the radial-velocity system is un- 
certain, there is very good agreement between the systemic (mass- 
center) velocity obtained separately from the individual datasets 
(OCA and CTIO data). 

Optical spectra of HD 164816 do not contain many features (see 
Fig. [5}: in addition to the interstellar lines there are only lines of 
the Balmer series and those of neutral and ionized helium and pos- 
sibly faint lines of CTL CIII, Si IV, and O II. In the OCA data for- 
bidden lines, e.g., [OIII] at 5006.84 A, are also visible, most likely 
residual features from the M8 nebular emission remaining after the 
background subtraction. The Balmer series lines (and also some He 
I lines) show narrow and variable nebula emission lines always at 
the same wavelength (e.g. 5876 A(see Fig. [6}). 

Our new spectra confirm the binarity (Penny 1996, Howarth et al. 
1997, Mason et al. 2009) of the system (see Fig.[6]l. We extracted 
the Doppler information and determined the spectroscopic orbij^] 
we used nine spectral regions centered at sufficiently strong He I/II 
lines (U 4024/4026, 4387, 4471, 4686, 4713, 4923, 5016, 5875, 
6678 A). The strong forbidden line [OIII] 5006.84 A has been re- 
moved prior to the fitting (visible mainly in the OCA data). 

The global modeling of the data included orbital parameters (P, 
e, a>, T , Vo, K 12 ), relative intensities an d projected rotational 
velocities of the components (vi^sini); see Table [4] for the indi- 
vidual parameters. Prior to the fit all emission seen in He lines has 
been cleaned "manually" in order to secure a proper determination 
of the radial-velocity. The fitted spectrum consisted of the sum of 
two appropriately broadened and Doppler shifted synthetic spec- 
tra for each orbital phase. Theoretical rotational profiles (see Gray, 
1976) were computed assuming a linear limb darkening law and 
solid-body rotation. For T e ff = 32500 K, logg = 4.00 [log(cm s" 2 )] 

3 http://www.astro.sunysb.edu/fwalter/SMARTS/ech_proc.txt 

4 In the case of the CTIO spectra only 4923 A, 5016 A, 5875 A, 6678 
A lines were covered 
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Figure 5. Average spectrum of all OCA observations on HD 164816 with 
line identifications (I.S. stands for interstellar band/line) 



and the wavelength range of the modeled He lines the linear limb 
darkening coefficient is 0.20 < u < 0.28 (see van Hamme, 1993). 
Because the shape of the rotational profile depends on the limb 
darkening only slightly, an average value of u = 0.24 was used. The 
template high-resolution synthetic spectrum taken from the Pollux 
database|^]corresponding to T e ff = 32500 K, log g = 4.00 and solar 
metallicity provided a very good fit to the data. Synthetic spectra 
corresponding to higher temperatures showed helium lines being 
partially in emission. On the other hand, a synthetic spectrum cor- 
responding to 30000 K resulted in smaller equivalent widths for 
the helium lines as observed. It is also possible that the helium 
abundance is higher than assumed in the synthetic spectra while 
the temperature is slightly lower. Previous determinations of the 
spectral type of HD164816 are mostly B0V (see e.g. Morris 1961). 
Our spectroscopy indicates 09V or 09.5V as the best estimates. 
In the first step we tried to find the approximate orbital period for 
the system using the full width at half maximum changes of the 
strong He I line at 5875 A. The present set of 29 spectra spanning 
17 months indicated only one possible orbital period, P ~ 3.82 
days. First fitting experiments, however, showed it was impossible 
to satisfactorily fit all the data without the assumption of a fast ap- 
sidal motion. Hence we added duj/dt to the parameter set with tj 
valid for the periastron passage Tq (at an arbitrarily selected epoch 
of the binary). 

Because the long-term All-Sky Automated Survey (ASAS^J light 
curve (Pojmanski, 1997 & 1998) does not show eclipses, a con- 
stant relative intensity of the components throughout orbital phases 
has been assumed. The best global fit to 29 spectra and 9/4 
(OCA/CTIO) spectral regions resulted in (among other parameters) 



5 http://pollux.graal.univ-montp 2.fr/ 

6 Please refer to http://www.astrouw.edu.pl/asas/ 
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Figure 6. The He I 5876 A line profile changes with the orbital revolution 
of the binary while the interstellar neutral sodium doublet (AA 5893, 5896A) 
is not following the orbital revolution. The same holds for the narrow nebula 
emission at 5876 A Hel, which is variable in strength but constant in wave- 
length, i.e. radial velocity. The dates of observations (in 2009) and orbital 
phases (counted from the periastron passage) are also shown. 



P = 3.81932(39) days and e = 0.232(16). The reduced X 2 = 1-394 
indicates a small inconsistence of the model and the fitted spec- 
tra. The larger^ 2 very probably arises from the simple assumption 
of spherical shapes for the components and also variations in the 
continuum rectification level. 

The resulting projected masses of the components mi sin 3 i = 
2.31(17) M Q and m 2 sin 3 1 = 2.10(16) M G indicate a low inclination 
angle for the system, which is supported by the lack of eclipses. 
The ZAMS mass of an 09.5V star (X = 0.70, Y = 0.28, Z = 0.02; 
see Claret 2004) is about 16 M Q ; the ZAMS radius is about 5 Rq (at 
the cluster age of 2.3 Myr (Prisinzano et al. 2005) the star would 
have a corresponding radius of about 5.6 R G ). Then the inclination 
angle would be 30 - 35°. Assuming i = 30°and a sin i = 16.8(2) R Q 
we get a semi-major axis of 33.6(4) R Q . In the case that the com- 
ponents are rotating synchronously (the synchronization time scale 
is usually two orders of magnitude shorter than circularization time 
scale; see e.g., Pan et al. 1998) their radii (using vifi sin i) are then 
R[ = 13.2 Rq and R 2 = 12.1 R Q , indicating that both components 
already left the main sequence or that they are still rotating asyn- 
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Table 4. Spectroscopic elements inferred from fitting the 29 OCA/CTIO 
spectra. The Icr error of the last digit is given in parentheses. Reduced^ 2 is 
given 



Parameter 



p 


[days] 


3.81932(39) 


e 




0.232(16) 


(OO 


[rad] 


-0.97(7) 


Vy 


[kms->] 


-77.1(13) 


A', 


[kms 4 ] 


109.4(26) 


K 2 


[kms^'l 


120.3(30) 


To 


[HJD] 


2455 000.88(3) 


dai/dt 


[rad/yr] 


0.67(7) 


U 


[yr] 


9.4(11) 


h 




0.536(14) 


V\ sin (' 


[kms^'l 


85.4(28) 


h 




0.463(14) 


V2 sin i 


[kms 4 ] 


79.9(32) 


mi sin 3 i 


[M G ] 


2.31(17) 


ni2 sin 3 i 


[M Q ] 


2.10(16) 


a sin i 


[Rq] 


16.8(2) 






1.394 



The following elements were optimized: P - period, e - eccen- 
tricity, <x>o - longitude of periastron passage at To, Vo - systemic 
velocity, K12 - semi-amplitudes of radial velocity changes, To 
- time of the periastron passage, dui/dt - apsidal motion rate (U 
is apsidal-motion period), I\g - relative intensities of spectra, 
vi ^sin/ - projected rotational velocities of components. Pro- 
jected masses of the components, m\£ sin 3 (', and projected ma- 
jor axis, a sin i are also given 

chronously (with rotation factors F l 2 ~ 2). Without knowledge of 
the inclination angle and in the view of possible asynchronous ro- 
tation of the components determination of the true radii and masses 
is also complicated. 

3.2 Distance 

The individual distance towards HD 164816 was determined by 
Megier et al. (2009) from interstellar Ca II absorption to be 864 ± 
136 pc, while the NGC 6530 cluster is at 1543 ± 345 pc (see Sect. 
1). Since these two values are deviant by more than 1 cr, we will 
discuss the distance of HD 164816 and whether it is or was a mem- 
ber of the NGC 6530 cluster. 

Since the radii of the two 09.5 stars in HD 164816 could not 
be determined directly by us, we cannot compute the distance by 
the Stefan-Boltzmann law from temperature, radius, and luminos- 
ity. However, since the two stars are located close to the dwarf se- 
quence (V), we can compute the distance modulus from the ap- 
parent magnitude of HD 164816 (corrected for both extinction and 
binarity) and the typical absolute magnitude of an 09.5V star. 

The published magnitudes of HD 164816 are (all in mag from 
Simbad and references therein) U = 6.22, B = 7.09, V = 7.09, I c 
= 6.99, J = 7.006, H = 8.053, and K = 7.072 (BV from Hipparcos, 
U from Reed et al. 2003, I Cousins from Rauw et al. 2002, JHK 
from 2MASS, Cutri et al. 2003), all having small errors of roughly 
±0.01 mag for UBVI and ±0.025 mag for JHK. 

Using the intrinsic UBVI C JHK colors of an 09.5V star ac- 
cording to Bessell et al. (1998) for temperatures of 31250 ± 150 
K and logg = 5.0 (for main sequence dwarfs), we can obtain the 
extinction by comparing the apparent and intrinsic colors; for the 
interstellar extinction law, we interpolate in each band according to 



Cardelli et al. 1989, Savage & Mathis 1979, and Rieke & Lebof- 
sky 1985. Then, we obtain A v = 1.10 ± 0.05 mag as interstellar 
extinction towards HD 164816. 

For the total bolometric luminosity of an 09.5V star, we use 
the latest determination from Hohle et al. (2010) using Hippar- 
cos distances and extinction corrections from Hipparcos BV and 
2MASS JHK colors, all corrected for multiplicity, interpolating be- 
tween 09V and BOV For the bolometric correction, we use again 
Bessell et al. (1998) for temperatures of 3 1250 ± 150Kandlogg = 
5.0 (for main sequence dwarfs), namely B.Cv = 3.025 ± 0.085 
mag. 

Then, we obtain as (main-sequence spectro-photometric) dis- 
tance towards HD 164816 the value 1030 ± 230 pc. This value is 
consistent with the value by Megier et al. (2009) from interstellar 
Ca II absorption being 864 ±136 pc, Hence, our assumptions ap- 
pear justified. For 09IV, 09V, B0IV and BOV stars, the values lie 
between 890 ± 40 pc and 1240 ± 140 pc, i.e. are consistent within 
the errors with the 09.5V case. In the case of the latter distance, 
the evidence for HD 164816 to lie in front of the cluster, would be 
much weaker. 

We conclude that HD 164816 may lie up to a few hundred 
pc in front of the NGC 6530 cluster. However, we cannot exclude 
that the star is located inside or at the front of the cluster. The pe- 
culiar radial velocity (corrected for Galactic rotation and solar mo- 
tion using a local standard of rest of (u, v, w) Q = (10.4, 11.6,6.1) 
km/s (Tetzlaffetal. 2011))ofHD 164816 is v rj ,„. = -80.7^°km/s 
whereas for NGC 6530 it is v rpec = -7.9^ ' km/s (Kharchenko et 
al. 2005). Hence, HD 164816 appears to move towards us relative 
to the NGC 6530 cluster. The difference between the two velocities 
is 72.8±7.9 km/s. While the radial velocity difference indicates that 
the star is now moving towards us relativ to the cluster, it can still 
be located inside or at the front edge of the cluster. If the system 
includes a neutron star born in a supernova, this supernova should 
have given the system a kick, which may have been the cause for the 
discrepant radial velocity. Given that the most massive and earliest 
star in this cluster (9 Sgr) has a spectral type of 04, the progenitor 
of the neutron star in HD 164816 has to have had an even earlier 
spectral type and, hence, had a life-time of below 3 Myr. Given 
that this is comparable to the cluster age (2.6 Myr), the presumable 
supernova should have happened very recently and the system HD 
164816 therefore cannot be located much foreground to the cluster. 
The proper motion of HD 164816 actually agrees with the typical 
proper motion of the NGC 6530 cluster. 



3.3 X-ray longterm variability 

As we have discovered a ~3.81 d orbital period we searched the 
available X-ray data sets mentioned in §2 for any longterm vari- 
ability; in particular those are three Chandra and one XMM-Newton 
set. 

For this purpose we first derived the absorbed fluxes for XMM- 
Newton and Chandra observations (see §2. 1 and §2.2) and plotted 
them vs the starting time of each observation taken from table [T] 
As there is only very limited statistics in each of the three Chandra 
exposures, we merge them and use the starting time of the middle 
observation. As a result one can see that the fluxes are consistent 
within the error bars. We have as well investigated the longterm 
behavior of the individual Chandra count rates, which cannot be 
compared to the XMM-Newton count rate as their instrumental re- 
sponses are different. We cannot find any variation in those count 
rates as well. Thus the merging of the datasets is justified. We hence 
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cannot detect the ~3.81 d orbital period in the available X-ray ob- 
servations. 



4 RADIO AND y-RAY OBSERVATIONS 

In order to get a complete view of the multiple system 
HD 164816 we have searched as well for any significant detec- 
tion in the available radio data catalogs. We could neither detect it 
in the National Radio Astronomy Observatory (NRAO) Very Large 
Array (VLA) Sky Survey (NVSS) nor by the GTEE 35 MHz Radio 
survey carried out by the low frequency T-array near Gauribidanur, 
India. Some feature could be detected in the H I All-Sky survey, 
the CO galactic plane survey and the 4850 MHz Survey carried 
out by the Parks and Green Bank Radio observatories. Due to the 
rather large pixel scales of around one degree, 0. 1 degree and 0.02 
degree respectively we consider the detected feature unlikely to be 
associated with HD 164816. In addition HD 164816 with a galactic 
latitude of — 1° is lying quite close to the galactic plane and the 
detected radio emission is hence more likely to come from diffuse 
emission of interstellar gas. 

With the FERMI y-ray telescope now in orbit we took as well 
data obtained by its Large Area Telescope (LAT) at the position 
of HD 164816 into consideration. As the whole sky is monitored 
every six hours for a bit more than two years now we analyzed 
data taken over this whole range (from 2008 Aug 04 to 2010 Aug 
09) in order to get the highest possible statistics. Events in the 
energy range from 600 MeV to 300 GeV have been chosen and 
have been filtered for galactic and extra galactic diffuse y-ray emis- 
sion. A likelihood analysis has been carried out afterwards which 
serves as a source detection since the photon statistics at the po- 
sition of HD 164816 are still very low. For this purpose spec- 
tral models for all sources in a ROI of ten degrees radius around 
HD 164816 have been accounted for contribution and were mod- 
eled simultaneously. A super exponential cutoff powerlaw has been 
assumed for HD 1 648 1 6 which is typical for y-ray emitting neutron 
stars (c.f. Trepl et al. 2010). Carrying out the likelihood analysis 
yielded no significant detection of the source. 



5 DISCUSSION 

The X-ray source associated with HD 164816 shows X-ray 
pulsations at ~ 9.78 s and in addition a soft X-ray excess with a 
blackbody temperature of ~ 49 eV both consistent with a compact 
companion (NS) in the system. Further optical observations 
confirm that HD 164816 is a spectroscopic binary consisting of 
two 09V or 09. 5V stars that are in a tight orbit of ~ 3.82 days. 
Accreting neutron stars and generally neutron stars in high-mass 
and supergiant X-ray binaries have in common that they emit at 
higher energies than 2 keV, but in the EPIC PN observation we did 
not detect any signal at above 2 keV, so that we can exclude an 
accreting pulsar as the source for the X-ray emission. 
All the aforementioned scenarios have in common that the NS is 
emitting non-thermal radiation that can be fitted by a powerlaw 
or Bremsstrahlung-model and is thought to occur from the ac- 
cretion process (either Roche-lobe overflow or wind accretion). 
Constant mass accretion can in fact be ruled out in the case of 
HD 164816 as the detected soft X-ray emission excess from 
2XMM J180356.8-241845 can be described by a pure blackbody 
model which means we observe thermal emission from the NS 
surface itself. In addition, evolutionary considerations based on 



standard formulae (Lipunov 1992), demonstrate that the stage of 
accretion is improbable taking into account the spin period of the 
NS and the parameters of the binary. Without better knowledge 
about the separation of the NS from the O-stars it is impossible to 
make clearer statements, but at a few Myr the NS is either in the 
ejector, or in the propeller stage. 

The orbital separation between the compact companion and 
the binary system must be at least around 100 Rq as empirical 
studies show that a tighter system is dynamically unstable unless 
the separation of the compact object from the binary center of 
masses is about three to five times larger than the semi-major axis 
of the binary system for an eccentricity of ~ 0.2 (Moriwaki & 
Nakagawa 2002) i.e. > 51 to 85 Rq here. Therefore wind accretion 
is negligible as the orbital separation is large. 
The other possible nature of the compact companion to 
HD 164816 might be a White Dwarf (WD) as the detected 
temperature of ~ 50 eV (5.80x 10 5 K) would place the object in the 
somewhat overlapping region of WDs as Super-Soft X-ray sources 
(SSS) in Be/X-ray binaries (Kahabka et al. 2006) and thermally 
emitting NS. 

However as the typical age of a WD should be > 1 Gyr the 
accompanying O-stars with masses of around 18 M Q and 20 M Q 
respectively would have already evolved and exploded. 
The compact companion lying just by chance in the line of sight 
of HD 164816 cannot be ruled out yet as the available Chandra 
observations are only spread by a few days (see Table [TJ. This 
is too short to significantly detect any separation between the 
centroid of the soft and hard component of the X-ray spectrum. 
If this detection is a chance projection, then the most probable 
candidates are radio pulsars. The probability for chance alignment 
of a neutron star within the point spread function (PSF) of the 
X-ray source can be estimated as follows: Neutron stars can be 
detected as either X-ray, y-ray or radio pulsars if they lie above the 
so-called dead-line in the P - P diagram; they reach this dead-line 
at an age of roughly 10 Myr (despite according to general cooling 
curves a NS may cool down to the temperature of ~50 eV at the 
age of ~5xl0 5 yr); if there are 5 x 10 s neutron stars in the Galaxy 
(at ages up to 12 Gyr), then there are roughly 3xl0 5 neutron stars 
detectable as either X-ray or radio pulsars (~ 2000 of them are 
known). Then, given the PSF of the Chandra source (full width at 
half maximum 3.70 ±1.60 arcsec) or the PSF of the XMM source 
(full width at half maximum 4.19 ± 1.78 arcsec) and the area of 
the whole sky, we obtain a probability for a chance alignment, 
between HD 164816 and a NS of 5 to 7 x 10~ 4 . However, the 
direction of these XMM and Chandra pointings is not random, but 
towards an OB cluster, where the probability for a NS is higher: 
According to Hohle et al. (2010), the area on the sky where almost 
all SNe are expected (inside OB associations), is 35 % of the total 
sky. The probability for chance alignment of a NS in the XMM and 
Chandra PSF is then 2 to 3 x 10~ 4 , still very low. Even if we restrict 
the estimate to the densest OB clusters close to the Galactic plane, 
the estimate would be less than one order of magnitude higher, so 
that the probability of chance alignment would still be below 1 %. 
Given that this estimate is very low, we have a large probability for 
the potential NS not to be a chance alignment, but to be related to 
HD 164816. This suggests that a NS is related to HD 164816. 
O stars are generally not known to have soft X-ray blackbody 
excesses. It has been shown by Naze 2009 that those stars are in 
general emitting spectra that can be fitted best with single or multi 
component MEKAL models with temperatures starting at ~0.2 
keV. 

A scenario that involves colliding winds instead of a compact com- 
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panion can be considered as well as HD 164816 consists of a close 
binary. However such a model would require a Raymond-Smith 
thermal plasma with a temperature of kT ~ 0.6 keV (Oskinova 
2005) or a hard non-thermal X-ray component (De Becker et al. 
2004). Taking a look at the MEKAL temperature of ~ 0.2 keV in 
the case of HD 164816 we note that it is almost a factor of three 
lower. In addition the spectrum is missing any component that can 
be fitted by a non-thermal model as the spectral counts from 2 keV 
on are consistent with zero. Hence, the NS does not appear to be 
accreting. 

As seen from §4 HD 164816 is not significant at radio wavelengths; 
this might have implications for the colliding wind scenario (De 
Becker et al. 2004), but that is beyond the scope of this paper. 
Assuming a distance 864± 136 pc we compute the X-ray luminosity 
to L x ~ 2.58xl0 31 erg s . Together with a bolometric luminosity 
for a 09.5V star from Hohle et al. (2010) this yields log(L x /L Bo/ ) = 
-6.73 which is in agreement with the relation for O stars found by 
Berghofer et al. (1997) (c.f. Fig. 4 therein). In the case of colliding 
winds this ratio should however be higher than the given value (see 
e.g. De Becker et al. 2004). It seems thus rather unlikely to get 
the low energy excess from colliding wind interaction. A similar 
conclusion was reached by Rauw et al. (2002), where they state 
that regarding the hxfL Bo i ratio no evidence for an increased value 
is found. 

Like Rauw et al. (2002) we find a higher than average interstellar 
absorption (see table |2]l which can be attributed to circumstellar 
or nebula material. This should result in strong suppression of 
the spectrum in the energy range 0.2 - 0.5 keV. In contrast to this 
we find an excess in that range on top of the absorption. Hence 
an additional source outside the region of circumstellar material 
responsible for that is highly probable. This is consistent with the 
earlier conclusion that the separation of the probable NS has to be 
at least 100 R . 

Emission line features in the X-ray spectrum as described in van 
der Meer et al. (2005) are not detected and as well not expected 
as such a phenomenon is only common in HMXBs and SGXBs. 
Those cases can be excluded as they both involve accretion 
processes which are not seen in our case. 

If there is a pulsating NS within the X-ray PSF and if this NS is 
not orbiting the O star of HD 164816, then it could be an isolated 
NS in the NGS 6530 cluster. A pulsation period of some 10s is 
not atypical for isolated young to intermediate-age NS - as e.g. the 
Magnificent Seven with pulsation periods of few to some 10s (see 
e.g. Haberl 2007). 



6 CONCLUSION 

We found peculiarities in the X-ray, optical, and kinematical data 
of HD 164816: 

• There is a soft excess in the X-ray spectrum of both XMM and 
Chandra. 

• The radius of the emitting (circular) area related to the soft 
excess would be ~7 km. 

• There is an indication of 9.78 s pulsation in the XMM data 
(not detectable with Chandra due to low timing resolution of the 
used ACIS-I detector (~3.24 s) and small count rate). 

• If HD 184816 includes a neutron star born in a supernova, this 
supernova should have given the system a kick, which is consistent 
with the fact that the star HD 164816 has a significantly different 
radial velocity than the cluster mean. 



All those four observational indications would be consistent 
with a compact object like a neutron star in the system orbiting 
the spectroscopic binary HD 164816 at a larger separation. A NS 
should have soft X-ray emission. A neutron star of few Myr age 
isolated from the other stars without accretion can have a period 
of few to 10 s (like the Magnificent Seven neutron stars). If there 
is a NS in the system, there should have been a supernova in the 
system before. Such a supernova could introduce a kick velocity to 
the neutron star born in the supernova. If the neutron star remained 
bound to the close spectroscopic binary, they would now all move 
away fast relative to their birth cluster. 

We did not find any evidence for X-ray emission due to collid- 
ing winds. An accreting neutron star would show hard X-ray emis- 
sion, which was not found. Hence, if there is a neutron star in the 
system it is effectively isolated (and non-accreting) from the O-type 
star. 
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